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Abstract 

A matrix model describing surfaces embedded in a Bethe lattice is considered. 
From the mean field point of view, it is equivalent to the Kazakov-Migdal 
induced gauge theory and therefore, ai N = oo and d > 1, the latter can 
be interpreted as a matrix model for infinite-tension strings. We show that, 
in the naive continuum limit, it is governed by the one-matrix-model saddle 
point with an upside-down potential. To derive mean field equations, we 
consider the one-matrix model in external field. As a simple application, its 
explicit solution in the case of the inverted W potential is given. 



1 Introduction 



Recently the "induced QCD" model proposed by V.Kazakov and A.Migdal [|l| 
has drawn much attention p|, |, H, |^. This model possesses the local ?7(1)®^ 
symmetry reminiscent to matrix models 0] which makes it inequivalent to 
QCD. Actually, to see that, it is sufficient to consider only the discrete Zjv 
subgroup which leads to the local- confinement selection rule incompatible 
with the usual QCD area law. 

However, as will be shown in the present paper, the Kazakov-Migdal 
(KM) model is interesting by itself. In the planar limit it is soluble by the 
mean field technique. But, within this approximation, it is equivalent to a 
matrix model with a Bethe-lattice (BL) embedding space. It means that all 
possible configurations are tree-like, i.e., have no embedding area. Therefore, 
it is natural to consider the KM model as an infinite-tension limit of a d- 
dimensional matrix model. Indeed, in this limit, surfaces degenerate to trees 
(or branched polymers in the other terms). Different aspects of branched 
polymer physics were discussed from the viewpoint of random surfaces in 
ref. 0. 

In some sense, this model is as interesting as any other non-critical string 
theory solved so far, because all of them do not have transverse modes on 
the world sheet. Moreover, except one-dimensional string theory, they do not 
have a real embedding space and can be interpreted only as the interaction 
of 2-d quantum gravity with different matter fields. But the d = 1 matrix 
model describes trees, because the world sheet has no area by construction. 
This model was solved in the planar limit in ref. ^ and in the double 
scaling in |jlO[ . The compactified embedding space was considered in ref. [p| 



where the exact solution for the singlet sector was obtained and interpreted 
as vortex-free string theory on a ring. Non-singlet sectors were considered in 
refs. [0, 0. 

Our mean field approach is a direct generalization of the one of refs. ||^, ^]. 
In the Hamiltonian formalism, free fermions interacting through an effective 
potential appear naturally and the Thomas-Fermi approximation is applica- 
ble. We consider only the planar limit, where the connection between KM 
and BL models holds and the mean field approximation is exact. In the 
continuum limit the Hamiltonian approach should be equivalent to the La- 
grangian one. The latter has been intensively explored in the induced QCD 
framework. Saddle point equations for the KM model were derived in ref. p| 
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and investigated in refs. [0, The exact solution for the quadratic potential 
was found in ref. 0. 

Mean field gives an exact answer usually for d bigger than some critical 
dimension and matrix models should not be an exception. It means that, in 
principle, we can use the BL matrix model as a convenient representation of 
multi-dimensional non-critical string theory. 

The outline of the present paper is the following. In Section 2, we establish 
the connection between KM and BL models. In Section 3, we describe the 
mean field approximation within the Hamiltonian formalism. To derive mean 
field equations, we consider in Section 4 the external field problem for a one- 
matrix integral and generalize the Brezin- Gross approach to it [|14|. So far 
only a solution for the cubic potential has been known |TB|. Our method 
allows us, in principle, to consider a general potential. In the planar limit, our 
approach gives a set of equations coinciding with the Migdal's one. However, 
a meaning of quantities involved is slightly different and, for finite A^, we 
still have rather simple linear equations convenient for the 1/N expansion. 
In Section 5, we investigate the continuum limit of the KM model. Some 
problems are discussed in Section 6. An explicit solution to the external field 
problem for the inverted W potential is given in Appendix. 



2 The KM model and Bethe lattices 

Let us start from a general matrix model with an embedding space being an 
arbitrary graph Q defined by its incidence matrix G^y {Gxy = 1, if there is a 
link connecting vertices x and y, and 0, otherwise). The matrix model action 
reads 

S = -NY^GxytrmxMy)) + Nj^tr (1) 

xy X 

where $(a;) is an A^ x A^ hermitian matrix attached to an a;'th vertex; the 
potential, Vo(A), is an arbitrary polynomial. The partition function is defined 
as the integral over all field configurations 
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In the planar limit (A^ oo), we adjust coefficients of Vo(A) so that to find 
a leading singularity of Z; universal continuum behavior takes place in its 
vicinity. At each vertex, can be decomposed into diagonal (f{x) and 

angular S{x) G U{N) parts: 

<^{x) = S+{x)^{x)S{x) (3) 
It is convenient to introduce the on-link gauge variables 

n,y = S+{x)S{y) (4) 
Obviously they obey the constraint that, for every closed loop L, 

n n,y = QL = i (5) 

{xy)eL 

where the product runs along L. 

If the graph ^ is a tree, all gauge variables are independent and can be 
integrated out by the Itzykson-Zuber (IZ) formula [p!6|: 



r J. dpf up^'t''^'^b 

li^, dn e^*^'^^'^^^ = iV-^(^-i)/2 TT n rf;\\^ , (6) 

n=l ^(W^W 

In eq. (j^), and ip without loss of generality are real and diagonal; A(0) = 
Y[i<j{'Pi ~ 0i) is the Van-der- Monde determinant. 

As a result, one gets a model where a role of dynamical variables is played 
by A^ eigenvalues of In the N ^ oo limit, the mean field approximation 

is applicable and, in this sense, the model is soluble. 

If the graph Q has loops, constraints (|^) can be imposed by averaging a 
number of 5-functions with the corresponding KM action: 



{L} {xy)&L 

where Z^^^^ is the KM partition function for Q. More precisely, we treat 
variables @) as ordinary gauge field, defined on a given lattice Q. To recover 
the original matrix model, this field has to be a pure gauge. However, eq. (j^) 
does not coincide with the weak coupling limit of a lattice gauge partition 
function, because the latter is singular and we have to impose the constraints 
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along non-contractible loops as well (if Q is not simply connected as a d- 
dimensional lattice). 

Expanding 5-functions in eq. (|^) in representations, we obtain the fol- 
lowing integral at each link 

/5(0,^) = j dU D^iU) e^'^^^^^^ (8) 

where D^fj{U) is a matrix element of a U{N) irrep R; lower indices run over 
a representation space Vr: a,b = 1, . . . , rf/j {d^ is the dimension of Vr); </> 
and ip are real and diagonal. 

The integral in eq. goes actually over the right/left coset 

f/(l)®^\f/(iV)/?7(l)®^ rather than over the U(N) Haar measure, since the 
action in eq. is invariant under left and right shifts by diagonal unitary 
matrices 

Uab ^ t/a^e^^"" (9) 

As was shown in ref. [|I3|, the symmetry gives rise to the selection 
rule for the integral (H): 

N 

jf,(0,V^) = O, unless ^ = (10) 

k=l 

where m^, k = 1, . . . , N, are highest-weight components of an irrep R (for 
U{N) they are unrestricted). Actually, as was noticed for the KM model in 
ref. [0, in order to obtain eq. ([T^), it is sufficient to consider only transfor- 
mations from the center of SU (N) . 

Inside representation spaces, Vr, we also have the selection rule 

/f^(0,^)=O, unless a, 6 e \4°^ (11) 

where V^°'* is the subspace of Vr spanned by all zero- weight vectors |13[ . Its 
dimension is equal to 

<-rn^x-«(o. (12) 

"''J k=l 
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where Xr{^^°') is the character of an irrep R. ^ only for irreps obeying 
eq. (0). When N ^ oo, 

d^R = Oi^) (13) 

It should be stressed that the selection rule ([TT|) cannot be derived from 
the Z]\f symmetry of the action in (||), because the center acts onto representa- 
tion spaces globally. The general approach to the calculation of integrals (M) 



was recently given in ref. [|17 



We can look at the (i-dimensional KM model also from another point 
of view. The absence of conditions on f2-variables is reminiscent to matrix 
models with tree-hke embedding spaces. We can consider a Bethe lattice, 
i.e., an infinite tree having the same coordination number for all vertices. 
This lattice is a covering of a regular lattice. We can obtain the KM model 
identifying eigenvalues ip{x) which cover the same vertices. The U{N) in- 
tegration includes the sum over all permutations and it makes eigenvalues 
indistinguishable. But it is sufficient to demand the equality of correspond- 
ing momenta of their distributions. These constraints can be interpreted as 
conditions on correlators in the BL model. Actually, we restrict only a small 
part of local degrees of freedom. It suggests that, in the N ^ oo limit, the 
KM model describes trees. Indeed, the mean field approximation is exact in 
this limit in both models and they are indistinguishable within it as far as 
their partition functions are concerned. 

This analogy holds also for correlators. Let us consider, for example, the 
following one 

tr(<l>(xi)<l'(x2))tr($(?/i)<I)(j/2)) 



( tr((^(a;i)l]L,^,^(^(a;2)l^L.,.Jtr((^(yi)fiL^^^^(^(y2)^^L,2,J ) (14) 

In the KM model we have to connect eigenvalues at different vertices by 
paths -L^i^2 multiply fi-matrices along them: 

^L.,.,= n (15) 

Owing to the selection rule (0), all appearing loops encircle zero area, be- 
cause the first non-trivial irrep obeying eq. (|T0|) is the adjoint. In the induced 
QCD framework, this phenomenon is known as the local confinement. If there 
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are links through which paths from different traces go, we can obtain a non- 
trivial answer also for more general loop configurations (an example is shown 
in Fig. 1). However, because of the second selection rule (|TT|) and eq. (P^), 
it has subleading order in A^. In general, averaging f2-matrices, we always 
lose half powers of A^, and a correlator has a correct order in N only if the 
corresponding index loops are pure backtrackings. Therefore, in the planar 
limit, the KM model is identical to the BL one. 

The ambiguity of the choice of paths can be interpreted as follows. As 
far as planar graphs are concerned, the KM model is defined actually on 
a covering of a regular lattice, where it is unambiguous and all possible 
configurations are tree-like. 

In next orders in A^, we lose this interpretation, as those two models are 
not equivalent anymore. The geometrical reason is clear. Already a torus has 
non-contractible loops which can wrap around plaquettes without creating 
an embedding space area. 

3 Mean field in the Hamiltonian approach 

Let us consider the Hamiltonian formulation of the KM model. It means 
that we treat one of dimensions separately as a time. In this respect, it is 
a direct generalization of the standard approach to the d = 1 matrix model 

We consider the partition function (|^) with the action (|l|). After the 
change of variables (H), (|) let us drop constraints (|^) and then integrate 
over ^-matrices in the t-direction only. If we introduce the continuous time 
keeping a lattice structure in "space directions", A^ eigenvalues at each site 
of the [d — l)-dimensional lattice become fermions in complete analogy with 
refs. P, ^. We have a fermionic ground state and, in the N ^ oo limit, 
the Thomas- Fermi approximation is exact. Each fermion, ipk, moves in an 
effective potential induced by its neighbours 

U{^k) = Vo{ipk) - 2{d - l)\ogI[ip,ip]{ipk)\p{i,)=piv) (16) 

where I[ip, i']{fk)\p{'4>)=p(ifi) means the IZ integral as a function of a fc'th eigen- 
value (fk when both densities coincide: pi^ip) = p(v^)- As we are looking for 
the ground state in the planar limit, densities of eigenvalues can be taken 
homogeneous in the space. Fermions at different sites interact only through 
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a collective field, since the IZ integral depends only on traces of powers of ip 
and ip separately. Hence, we have no Fermi sphere and can consider fermions 
at each site independently. The one-particle Hamiltonian reads 

H{p,x) = ]^p' + U{x) (17) 

p and X are a momentum and a coordinate, respectively. Now, we can repeat 
standard steps [^. The Fermi level, Ep, is defined by the equation 

N = J^eiEp-H{p,x)) (18) 

where 9{x) is the step function. 

The ground state energy is given by the formula 

F = H{p, x) d{Ep - H{p, x)) (19) 

J Zn 

Differentiating eq. (0) with respect to Ep, we find the inverse frequency 
of a particle moving in the effective potential 

dN 1 f dpdx 



J^6{Ep-H{p,x)) (20) 



dEp Lj{Ep) 

A simple calculation gives the local density of eigenvalues 

r'^F dE 1 

pW^ ■-» ;;v^ ^_i_^/^^;r^ (21) 

Critical behavior occurs when the Fermi level reaches a local maximum 
of the effective potential U{x), which is, in general, quadratic: 

U{x) = U{xo)-^ix-Xoy + ... (22) 

However, in our case, there is a possibility to get an answer different from 
the d = 1 model, as the coefficient a can be a function of the Fermi level Ep. 
Let us introduce the analytic function 

fix) = f dyP^ (23) 
J X — y 
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where the integral goes over a support of p{y)- From the one- matrix-model so- 
lution we know that, in order to get a given density of eigenvalues, ^Im /(x), 
we have to take 2 /Re/(x) as the one- matrix potential. The IZ interactions 
between lattice sites perturb this simple solution. 

To obtain a closed set of equations, we have to express U{x) through p{x). 
One of possible ways to do it is to consider the one-matrix model in external 
field. 



4 External field problem revised 

Let us consider the matrix integral with the linear source term: 

I[X] = J rf^Vexp Ntr{XY - V{Y)} (24) 

where X and Y are N x N hermitian matrices, and V{Y) is an arbitrary 
function defined by its Taylor expansion V{y) = J2 VkV^. 

It is obvious that 1[X] depends only on eigenvalues of X. The standard 
way to deal with the integral (0) is to write down the Schwinger-Dyson 
equation 

j d''^Yii[x - V'{Y)\ exp Ntr{XY - V{Y)} = (25) 

and, then, rewrite it in terms of eigenvalues of X. 
It is convenient to introduce the resolvent 

Gij{z) = J ^^).^.exp NtT{XY - V(Y)} (26) 

Then, eq. (^) can be rewritten as follows 

dz 



■tr 
2m 



{X -V'{z))G{z) =0 (27) 



The integral goes along a small circle around z = 0. 

Eq. (pTl) has to be accompanied by an equation for G{z) which can be 
obtained from the following obvious identity 

X|A1 = / l^)(-^)}e-<x.-v-,.„ (28) 
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Let us take X almost diagonal 

X = x + i[A,x\ + ... (29) 
where A is an infinitesimal hermitian matrix: ||y4|| ^ 1. [ , ] is the commu- 



tator. Then, eq. (|28| ) can be rewritten in the form 



1 r 

d-Y (( 



1 5 11 

z - — T^ldij + 



1 

Z — Y '^3 

1 9 



d 



N Xj — Xj id A. 



z-Y 



^Ntr{xY-V{Y)) 



A=0 



Ndxi 



where 



Ntr{xY-V(Y)) 



(30) 



(31) 



are diagonal elements of the resolvent matrix Gij{z). If we normalize them 

as 



then we find the system of equations 



9k[z) 
X[X] 



dz 

Xk=f —V'iz)Wkiz) 



(32) 



(33) 



Eq. (|3^ gives the following Laurent expansion 
and, hence, we can recover T[X] from a residue of zWkiz) 



(35) 
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Eqs. (^), (0) are convenient for the large limit to be taken and also 
provide a suitable framework for the 1/N expansion. 

Let us suppose that, in the ^ oo limit, there exists the density of 
eigenvalues having a finite support 

1 ^ 

p(^) = Ji5;L^E^(3^-3;fc) (36) 
fc=l 



Following ref. we can introduce a continuous variable x instead of the 
lower index of the function Wfc(-2): 'w{z,x). Then, in the iV — >■ oo limit, eqs. 
31, (H take the form 

x = j ^y\z)w{z,x) (37) 

1 = {z- Wi{x))w{z, x)~ dyp[y) (38) 

J X y 



where 



is the second coefficient in the Laurent expansion 

^(^'^) = E-Fi:r (40) 

k=0 ^ 

It should be noted that eq. (|38| ) coincides with the one obtained in ref. in 
a slightly different context. 

These equations are valid, strictly speaking, only on a support of p{x). 
Let us suppose, following ref. Ijl^, that they can be regarded as dispersion 
relations and, in such a way, can be given a meaning on the whole complex 
plane. 

Let us introduce the function 

f(,,,)=/,,^toMM) = f AW (41) 
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where the integral runs over a support of p{x). F{z,x) is analytic at the 
infinity and 



Im F{z, x) = —p{x)w{z,x) (42) 

TT 

In terms of residues, eq. ( |5B| ) can be written as follows 

Wk+i{x) = (wi(x) + fo{x))wk{x) - fk{x) (43) 

The functions fk{x) are not known (except fo{x), which is determined by 
p{x)). However, if the potential V{z) is a polynomial, we have a finite number 
of equations and a bootstrap solution can be, in principle, found, provided 
an analytic structure of all functions is guessed. No further information is 
needed. In this respect, our method is a direct generalization of the Brezin- 
Gross technique P |. 



A singularity of Wi{x) at the infinity is determined by the equation 
V'{wi{x)) = X. Therefore, for potentials of a degree bigger than 5, this 
singularity is not algebraic and the problem is very complicated. However, 
for the inverted W potential, wi{x) can be found in a closed form by using 
the well-known solution to the cubic equation. For reader's convenience we 
give it in Appendix. 



5 Continuum limit and scaling at d > 1 



In the case under consideration, the density of eigenvalues of external field 
coincides with the one of the integrand and we have the following equality 



dx p{x) w{z, x) 



dy 



pjy) 

z-y 



(44) 



In ref. Migdal solved eq. ( pHD for w{z^x) after expressing wi{x) from 
the Lagrangian saddle-point equation 



Wi[X 



^Xix)--Refix) 



Then, the condition 
Re/(x; 



2d ' d 
gave the following Master Field Equation 



J TT 



Im log 



X 



-W'{y)+t\mf{y) 



(45) 



(46) 
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where the integral goes over a support of Im/(y); W'{y) is the derivative of 
the effective potential 



U{x) = Vo{x) - 2{d -1) J dx wi{x) (47) 
which, after the substitution (|^) for wi{x), takes the form 

W{x) = ^Vo{x) + ^^^^Y^ J ^^^(^) 

In the Hamiltonian approach we easily obtain the square root form of the 
density (pi]), while in the Lagrangian one the simple saddle-point equation 



(45) holds. Both approaches are compatible in the continuum limit. In the 
simplest, d=l, case it was demonstrated in ref. 0. Let us use the same 
technique also ior d > 1. 

To take the continuum limit, we have to rescale the field and the potential 

X = Xe'''^-^- Vo{x) = X^de''-^ + e'' v{X) (49) 

where £ — > is a lattice spacing. To preserve the normalization, we define 

f{x) = e'-'/^iX) (50) 
In terms of new variables, eq. (E6[) takes the form 



Re^(A) 



rda 



(51) 

In the £ ^ (or d oo) limit, the term e'^v'{fi) seems to be negligible 
with respect to e'^~'^Iieip{fi). However, for the sensible continuum limit to 
exist, Reip{X) should become 0{e'^) dynamically. 

Let us expand the r.h.s. of eq. ( PD in powers of e, 



Re^(A)=/^^^(l 

J IT A — /i I 



+ 



1 



X — fi \2d 
[ImV^(/^)]^ 

3{x-^^y 



d- 1 



+ ... 



(52) 



as well as v^(A) 
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Re(/^(A) = £"iri(A) + £"V2(A) + . . . 



lm^{\)=e^^J^{\) + ... (53) 
Substituting these expansions in eq. (|52D, we find 

v'{X) 



2(d- 1) 



K2 = d + 2 j^(X) = J^i^-ll I dXr^iX) (54) 



in obvious agreement with eq. (|21|). 

The leading term gives the equation for r2(A) 



2{d-l) J n A-/i 

It has the form of the one-matrix-model saddle-point equation with the 
upside-down potential. It is tempting to connect such a strange phenomenon 
with the appearance of a tachyon in multidimensional string theories. Its 
technical reason is clear: the effective potential tunes itself to cancel a bare 
one up to higher order terms. Though this branch corresponds formally to 
the continuum limit, it seems to describe an effective zero- dimensional sys- 
tem having pathological properties. It is unstable for most potentials, when 
the initial lattice model is obviously well-defined. 

If we keep only a quadratic term in the potential, v{fj,) = ^/x^, the 
solution is known from ref. to be 




4 - (56) 



In the e — > limit, there are two possible values of 



c 



rl 1 

em 
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and 



- 2d-l ' ^^^^ 

The latter, C2, does not correspond to the continuum hmit. However, this 
branch is always stable. If eq. ( |55[ ) has no solution, the model should become 
gaussian in the e — limit regardless of a form of the bare potential. Of 
course, one could interpret this saddle-point as a non-trivial background and 
construct an e expansion around it. But it is doubtful that this would have 
connection with continuum theory. This solution corresponds to the strong 
coupling phase of the lattice model. In the case of a quadratic potential it 
was discussed by Gross |Q. 

However, there is some hope to obtain non-trivial behavior at least in the 
Hamiltonian approach for d ^ 1. A possible critical scaling could correspond 
here to the "upside-down" effective potential (|23) in which case 



ImcfiX) ~ V-e + aA2 (59) 
In the scaling limit, e ^ 0, we have 

Reip(\) = 2 4^ V-g + a/i^ ^ -\^\og ( —) + 0(1) (60) 

J r-r A — U ^ e ^ 

y/e/a ' 

where we have introduced the cutoff A = 2A^/a. 

In the d ^ 00 limit, the presence of the continuous time should not 
influence an effective potential which is proportional (as follows from eq. (^8])) 
to 2 /Re/(x). If we assume that the quadratic top (^2]) is really created 
by non-linearities at large distances, then, from the form of the effective 
potential, we find the following selfconsistency equation 



from which it follows that 



up to log-log terms. 



A 2 

a = 2Va£'^log(— ) (61) 



2e^log(-))' (62) 
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Now, we are in a position to calculate the singularity of ^j^^jy- From 



eq. p^) we have e ~ U{xo) — Ep and, hence. 



d 


\ 1 1 


d f 


de 







' e/ a 



Finally, we find that 



dx . ~ — - (63) 

V -e + aa;2 e|loge| 



log I log e I (64) 



uj{Ep 

Repeating standard steps p, we can introduce a couphng constant, by 
rescaling x and find 

dF _ Ep _^ dg_ 

dEp ~ u{Ep)~'''dEp ^^^^ 

A critical value, gc, corresponds to Ep touching the top of the potential. 
Then, we get for the second derivative of the free energy with respect to the 
coupling constant the following formula 

^" ~ 1 11^ A I (66) 
log I log 5^ I 

where Sg = gc — g. Let us remind that we have assumed that d ^ 1. 

Unfortunately, we have no real derivation of this result, and it should 
be regarded only as a reasonable proposal. If this possibility fails, it will 
presumably mean that, in the continuum limit, the KM model is effectively 
gaussian (if we do not take eq. ( |53| ) seriously). 

For usual branched polymers, the second derivative of the free energy 
has a square root singularity, which corresponds to the string susceptibility 
exponent •jstr = 1/2 IQ- In our case •jstr = 0. The reason of it could be 
that random surfaces have internal degrees of freedom which produce some 
weights at branching points of polymers. The double logarithm should not 
be surprising, as for random surfaces •jstj- < 0, and, at ci = 1, we have 
already had F" ~ 1/ log6g. 
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6 Discussion 



The Bethe lattices are considered in statistical mechanics as simple systems 
for which mean field gives exact answers. In a sense, they correspond to 
the limit d oo, i.e., their effective dimension is infinite. It is known that, 
for local 0^-type interactions (for example), this regime is exact actually for 
all dimensions bigger than 4. Presumably, it should be the case for matrix 
models as well. If it is really so, the obtained mean field solution is exact for 
all dimensions above a critical one and large dimensional string theory indeed 
describes branched polymers as was argued in ref. 0. For the KM model 
itself, the critical dimension is obviously equal to 1. If it is independent of a 
tension, non-trivial string theory does not exist. 

However, the matrix models do not contain a bare string tension as a 
parameter. If we substitute the heat kernel 

Kp{Q)=j:dRe-f'^'^Xam (67) 

R 

for 5-functions in eq. {Cr is a second Casimir), we obtain a model in- 
terpolating between the matrix and KM models[|. Though we have no clear 
interpretation for intermediate values of (3, it may bear some properties of a 
tension dependent matrix model. Of course, its dynamics is extremely com- 
plicated, as we have a hybrid of a matrix model and lattice QCD. But, if 
QCD string really exists, it could be used to introduce a bare string tension 
in a matrix model. And, from this point of view, the interpretation of f3 as 
a tension may be quite reasonable. 

To conclude, we should say that the mean field problem for matrix models 
is still far from being well understood. Our attempt to construct a non-trivial 
continuum limit for the KM model at > 1 has obviously failed. At least, 
eq. (|55D does not seem to be good for it. On the other hand, eq. (0), though 
looks rather natural, was not rigorously derived. It is even not clear whether 
it really corresponds to a physical solution of the model. However, in order 
to prove (or disapprove) it, we would have to solve eq. (^) with a non-trivial 
potential, which does not look to be a simple task. Also an interpretation of 
the KM model in subleading orders in remains to be done. Even for the 
d = 1 model compactified on a circle, a clear understanding of non-singlet 

^Another proposal can be found in rcf. 
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sectors is still absent. Let us hope that subsequent works will clarify all these 
problems. 
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A Appendix 

We consider eq. (|2^) with V{y) substituted by 



V{y) = - -/ 



(68) 



Then eq. ( pT] ) takes the form 



A solution to this equation can be written in the form 
where 



(70) 



r X 



2 V ^3' 



and u{x) and v{x) have the following asymptotics 



u{x) = 1 + 0(l/x) 



v{x) = - + 0(l/x) 



(71) 



(72) 



From the condition that Imwi(x) = 0, when Im/o(x) 7^ 0, we find the only 
possible real solution 



Wi{x) = r^/^(x) + 



3^1/3^2;) 3 J X — y 



r{x)' 



-2 



(73) 



By substituting this ansatz in eq. (p9|), one can check after a tedious algebra 
that the former obeys the latter only on a support of p{x). For a cubic 
potential, the corresponding equation was satisfied on the whole complex 
plane |l^. The other peculiarity of the quartic solution is that we have no 
additional equation for a position of the cut. 

Integrating over x we find the formula for the integral (|^) in the planar 
limit 



N 



X + 



19 



log 



r X \ 3 



^(r^\x) + (^)V ^^\x)) + \J dxp{x) J dyp{y) 



r{y) 



+ 



r{x) 



+ 1 



r[x)\3 



r{y) 



+ 



r{y)Y 

r{x) 



+ 



+ log 
3 



^(r{x)r{y)y^^ + ^(r{x)r{y)) 



+ 1 



^(r(x)r(y))^/%^(r(x)r(y))-^/']} (74) 
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